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ABSTRACT 

A brief review of advances in the inverse design and optimization theory in the 
following engineering fields in China is presented: I) Turbomachine aerodynamic inverse 
design: including mainly: (1) two original approaches — image-space approach and varia- 
tional approach , (2) improved mean-streamline (stream surface) method, (3) optimization 
theory based on optimal control. II) Other engineering fields: inverse problem of heat 
conduction, free-surface flow, variational cogeneration of optimal grid and flow field, 
optimal meshing theory of gears. 


I. INTRODUCTION 

Up to now, mostCover 95*) of the technical literature deals only with the direct 
(analysis) problem due to possibly the fact that the inverse(desi gn) problem( finding the 
unknown boundary shape) is, in general, much more difficult to formulate as well as 
to solve than the direct one, though the inverse problem is more important for , and di 
rectly related to, practical design. As a result, for instance, almost all turbomachine 
bladings are still designed by repeated use of direct problem methods in a cut-and try 
manner, which is of course not only inconvenient and time-consuming, but also incapable 
of providing very good results. So in the 1950's in China we have tried to apply the mean- 
streamline method for inverse problem of Wu & Brown[30] to cascade design and some impro 
vements of this method were suggested [2] . In the 1960's a new image-plane approach to the 
inverse problem was proposed [5]. It was realized, however, that also the inverse problem 
can not be successfully used for practical blade design, because it often leads to blade 
configurations that are either unfeasible from consideration of stress, vibration, coo- 
ling and technology or even unreal izab le( e. g. giving profiles unclosed or with ne- 
gative thickness). Therefore the traditional direct and inverse problems can not keep 
up with the development of modern turbomachinery(TM) and it was suggested in Refs. [18, 56] 
to extend the scope of aerodynamic problems and reclassify them into four categories: di- 
rect, inverse, hybrid and optimization problems. Then the image-plane approach was ex- 
tended to hybrid problem in Refs. [6, 56], and another new approach to inverse and hy- 
brid problems based on variational principles (VPs) was also suggested in Refs. [18,21]. 
Since then, a lot of variants of the image-plane approach, the variational approach and 
the mean-streamline method have been developed in China and extended to 3-D case. 

The Chinese research on the optimization problem of bindings started with the pro- 
blem of optimal radial distribution of flow parameters in TM with long twisted blades in 
1963[52]. Later, advances in this area are characterized and facilitated considerably by 
the introduction of modern optimal control theory. 
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II. RECLASSIFICATION OF ENGINEERING PROBLEM SETTING 

Generally speaking, any problem of engineering sciences can be posed in different 
ways, resulting in four problem categories: direct, inverse, hybrid and optimization pro- 
blems. Specifically, for the aerodynamic problem of blade cascades these problems are 
defined conceptually in Table I. The aerodynamic problem for S 2 “Stream surface can be 
classified similarly as shown in Table II. 

The hybrid problem is a unification as well as a generalization of the direct and 
inverse problems, encompasses a wide variety of typesCsee Table III for cascades on 
arbitrary streamsheet of revolution) and hence is very flexible and capable of meeting 
various design requirements. It provides design engineers with a series of new rational 
versatile ways for blade design. In addition, the inverse and hybrid problems also cons- 
titute an important ingredient of the optimization problem. 


■ 


To be sought 


Table I. Problem Classification of Cascade Flow 


problem Given 


Direct (D) 


Inverse(I) surface velocity distribution cascade geometry 


cascade geometry 

velocity field 

________ 



3 Hybrid (H) 'partly geometric conditions' A the remainding unknown 
'partly flow condition' geometry * flow field 


4 Optimiza- Dobjective functional optimal cascade geometry 

tion (Opt) 2)design inequality constraints * flow field 


Table II. Problem Classification of S 2 -Flow 



Hybrid 


Optimization 


on hub * casing 


geometry 


pressure distribution 


pressure distribution 


geometry 


Shape or V y r 'partly geometry'+ 'part ly pressure 
distribution' 
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Table III. Hybrid Problem Types of Cascades 


Types 

Given condi tions 

geometric 

aerodynamic 


part of airfoil form 

velocity distribution along the 
remainding part of profile contour 

h b 

airfoil thickness 
distribution 

blade loading distribution 

<p p -p s ) 

H C 

airfoil thickness 
distribution 

velocity distribution along the 
suction surface 


airfoil thickness 
distribution 

distribution of velocity difference 
<W S -W p ) 


Since the hybrid problem of fully 3-D flow may have a wide variety of types* 
depending on the manner in which the boundary conditions (B.C.) on the blade surface are 
combined with those on the annuclar walls/ it is necessary to employ some properly 
defined compound symbols to designate them as proposed in Refs. [46/47], For instance* the 
symbol (IXH^) designates such a hybrid problem type in which an inverse problem is posed 
on the blade surface* while a H^-problem is posed on the annular wal ls(Fi 9 . la). In other 
words# the symbol before 'X' characterizes the problem type on the blade surface* while 

the symbol behind 'X' that on the annular walls. It is easy to see that the scope of 

possible hybrid problem types can be made even much broader by posing different problem 
types on different portions of the bladeCand/or annular) walls (Fig. lb). 

III. INVERSE & HYBRID PROBLEMS OF BLADE-TO-BLADE FLOW IN TURBOMACHINES 

The research on inverse and hybrid problems In China has been going basically along 
the following three lines: (i) universal approach based on image-plane concept; (ii) uni- 
fied approach based on variational principles (VPs) and the related finite element method 
(FEM); (iii) improvement of the mean-streamline method and of other well-known methods. 

1) Universal Approach Based on Image-Plane Concept 

Two different image planes ^ and %rj have been introduced: 

1-1. Methods based on image-plane . 

The first universal image-plane method for solving the inverse problem of 2 -D com- 
pressible cascade flow was suggested by Liu in 1964 [5] and extended to cascades on 
arbitrary streamsheet of revolution by Liu * Tao in 1967 [ 6 ] and to hybrid problem H 4 by 
Liu * Tao in 1981 [ 6 ]. The main difficulty is the treatment of unknown boundary(b lade 
surface) and was successfully overcome by introducing a nonor thogonal curvi linear(stream- 
line) coordinate system (von-Mises coordinates) defined by (Fig. 2> where *] should be rep- 
laced by Y ) 


ior-. 5=M)] 

r= r(i<f) 


where the stream function Y * s defined by 
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'dip' _ . ?)W 

9T = _t P/V> ~d<f = P^l (2) 

It is expedient to regard Eq.(l) as a Mapping* which transforms the original irregu- 
lar periodic flow domain with soae unknown boundaryCAB in H^-problem* Fig. 2a) on the 
physical stream surface into a simple rectangular one with fully known boundary in the 
iaage plane (Fig. 2b). Moreover* the following four alternative formulations were derived 
in Ref. [5* 6] : 

i) first-order partial differential equation (PDE) system 

i i) second-order PDE 

iii) integro~di f ferent ial equation 

iv) integral equation system* 

of which only the integro-differential formulation for homentropic flow is given here for 
reference [6]: 





( 5 ) 




< M) 


Here /(— /Cl 0 ) and /J u (= /Cl,) >re diaensionless relative velocity and blade speed 

respectively) m = (J C-l) 1 * a 0 reference speed of sound; cO — rotor angular speed; SC — 
the contour abeda of a finite area SA (Fig. 2b). 

A H^-and an inverse problems of a cascade on a general streaasurface of revolution 
have been solved by this aethod by Chen et al. in Ref. (7]. 

This aethod was then iaproved considerably by Liu [8] via introducing a new aoaent 
function fl defined by 


!f= Xr ^ , !$— (V/Jr <1> 

The aoaent function has soae special features, for instance! (i) its increaent 
around any closed contour enclosing an airfoil is just equal to the aerodynaaic ao~ 
aent M z exerted on airfoilt 


Afi= m 2 = r a f 


This is just a generalization of the wel l-known Kutta' s lift theorea for 2-D flow, 
showing that is proportional to both absolute circulation around airfoil f and flow 
rate through an interblade channel ATj* (ii) the pressure p can be coaputed directly froa 
Eq.(5>, resulting in two advantages: first, no density aabiguityllO] appears; second, for 
the inverse problea the B.C.Ci.e. the distribution of fl ) on the airfoil contour is of 
the Dirichtet's type and hence easy to deal with. Also in this case the four alternative 
foraulations aentioned above are possible, of which only the second-order PDE foraulation 
is given below: 
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( 8 ) 


where A/B/ C/D/E are functions of Mach number/ P / % and (Fig. 2); F is a function of the 
gradients of entropy and rothalpy. In Ref. [8] an inverse problem of a cascade on a 
conical s treaasur face taken froa Ref. [78] was solved and the result is given in Fig. 3. 

Later/ a nuaber of aethods using this iaage plane have been also published in 
Re f s [9—1 2] ^ differing froa one another/ however/ by different choice between the four 
above-aenti oned foraulations and by different iterative strategies. Thus/ in contrast to 
Refs. [5/6]/ Shen & Ma [9] solved the -and H c -problems of 2-D transonic cascade flow 
by eaploying the lst-order PDE foraulation and Jameson's rotated difference scheae/ while 
Chen * Zhang[10]/ using the second-order PDE foraulation for the dependent vari ab le 
presented a nuaerical method for solving direct/ inverse and hybrid -probleas along 
with three nuaerical inverse problem examples/ of which the one for a tandem cascade is 
given in Fig. 4. Some difference between the calculated and original profiles might be at- 
tributed to the use of the measured velocity distribution as input for the calculation. 
This aethod has been modified by Sun et al. in Ref. [11] by using a boundary-fitted coor- 
dinate (see Eq.C9>) instead of . The nuaerical result of a supercritical cascade to- 
gether with its modified design is shown in Figs. 5 * 6. In Ref. [12] a method similar to 
Ref. [10] for ~ * I-probleas was presented for rotational flow/ and a method for remov- 
ing the density ambiguity is also given. In addition/ a rational cascade design proce- 
dure consisting of successive use of H^, -and I-probleas is proposed. 

1-2. Methods based on inage-plane §*7 . 

All aethods using inage-plane suffer froa the shortcoming that singularities 

appear in the vicinity of blunt leading and trailing edges due to local aul tivaluedness 
of the napping Eq.(l). To circumvent this difficulty/ another aethod for hybrid problems 
was suggested by Liu [13]/ where a new inage-plane defined by (Fig. 2) 


1=1 Cor-. 5=f«)J 

i= vf 


was introduced/ where 7? and 7s < the T| -values on suction * pressure sides) are given 
constants. Also in this case four alternative foraulations can be derived/ but only 
the integro-di fferential foraulation is given here for reference (Fig. 2b)i 



f|=xpOVV''V)> ff=Ti°V.< 


(10) 


( 11 ) 

( 12 ) 


where tg'jf is the slope of the ^-coordinate linei H^= ^ ^A T }jT r }s) is the scale factor of 
the coordinate 1 | ; is the angular width of the blade channel. We can see that 
this new iaage-plane aethod is particularly advantageous for solving those hybrid prob- 
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leas with given airfoil thickness (and hence is also known). 

Other aethods based on i«age plane. using/ however, the following second-order 
PDE foraulation. have been presented by Chen et al. [14.15] and Ge [1 6] t 


A ./ ^ V_ + A tfV+A 22E+A QV— b 






where A,. A^. A 3 . A^. A s are functions of p and the aetric tensor g 2J . g J3 ; Bj depends 
gradients of entropy and rothalpy. velocity.^). (T and viscous forces. In Ref. [14] 
problea of potential flow is solved and one of the numerical examples is given in 
The figures 8 & 9 taken froa the viscous inverse problea solutions of Refs. [15] 
respectively show that for the saae inlet and outlet flow angles the airfoil in 
flow is aore strongly curved than that in inviscid flow. 


( 13) 

on the 
the H/|- 
Fi g. 7. 

* 116] 
viscous 


2) Unified Approach Based on Variational Principles 

Basically there have been developed two coapletely different variational approaches, 
following a systeaatic way suggested by Liu [17] . 

2-1. Approach based on VPs in the iaage plane 3lff(Fig.2) 

In Ref. [18] Liu established two faailies of VPs and generalized VPs in terns of the 
aoaent function XL and angular function respectively for the and H B -probtens in the 
iaage plane 5Y- which were nodified by Liu & Yao to give the VPs for the H^-problea in 
Ref. [19]. Only one of these VPs is given below: 

(£ff X" +hJ}<*iW+L Jl <H> 

where the boundary integral tern L^takes different for® for different problem types. 
Based on these VPs involving^/ some finite element (FE) solutions to H^-and H B -problems 
have been presented in Ref. [20] by Yao et al.* fro* which Figs. 10 £ 11 for a cascade on a 
conical stream surface [78] are taken. 

2~2. Approach based on VPs with variable domain. 

Making use of the functional variation with variable domain* Liu was able to 
establish three families of VPs and generalized VPs for Ify-* H B ~and H^-problems in terns 
of the potential and streaa functions § & t^/ for potential and rotational flows in Refs. 
[21*22] and extended thea to transonic flow with shocks in Ref. [24]. Moreover* variable- 
doaain VPs using Clebsch variables have been also developed for 2-D transonic rotational 
channel flow by Liu [25]. 

Numerical solutions to H^-and H^-probleas based on VPs of Refs. [21*22] have been 
obtained by Yan & Liu[22* 23] by aeans of a new finite eleaent with self-adjusting nodes 
for numerical realization of the functional variation with variable doaainCFigs. 1 2 * 13). 

Perhaps a very attractive merit of this variable-domain approach is that it can be 
straightforwardly extended to fully 3-D flow. 


3) The Mean-Streamline Method (MSLM). 

This method originally suggested by Wu * Brown [30] was improved in many aspects in 
China. A survey of this development before 1984 has been presented by Cai[2]. Recent 
research includes Cai's paper [3] and Wang's paper [4]. 


4) Miscellaneous Approaches. 

Several known approaches to inverse design of cascades were improved or modified in 
China. 

4-1. Iterative method based on direct problem solver. 

Such a method is suggested by Wang in Ref. [26] to solve inverse and various hybrid 
probtemsC including * H c ) and extended to viscous flow in Ref. [27] by incorporating a 
boundary layer solver of integral type. Based on this method* Wang et al. proposed a 
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quasi-3-D design procedure for impe l lers [28; 291 . 

4~2. Time-dependent Method. 

Starting from the integral fori of aerodynamic equations; a f ini te-volume method for 
inverse cascade problem is given by Zhou & Zhu [31]. 

4-3. Hodograph method. 

It was improved in the transonic region by incorporating some analytical nozzle so- 
lutions and generalized to cascade flow along general streamsheet of revolution indepen- 
dently by Chen [32] and Yao[33], 

IV. INVERSE & HYBRID PROBLEMS OF S 2 -FL0W IN TURBOMACHINES 

Similarly to Sj-flow; the image-plane approach and the VP-based approach mentioned 
above can be applied to S 2 ~flow as well. 

1) Unified VP-Based Approach. 

Starting from the basic equations of Wu ; s S 2 -flow mode l [T 79] » first complete VPs 
and generalized VPs for the semi-inverse problem were established by Liu [34] and the 
corresponding FE solutions were obtained by Qin et al. [35], Inverse and hybrid problems 
of S 2 -flow were formulated in a unified manner by VPs with variable domain by Liu [36] and 
by VPs in an image-plane ^ by Cai * Liu [37] ; which have been generalized to flow of 
pure substance by Xu [38], In Ref. [39] VPs for hybrid problems of axisymmetric channel 
flow were derived by Tao & Liu. 

2) Universial Image-Plane Approach. 

Using an image plane ant * 9iven a distribution of circulation V^f on S 2 ~surface> 

Ge presented a method for solving the complete inverse problem and a hybrid problemCwith 
unknown hubCor casing) wall* see Table II); thereby a second-order PDE for was 

derived and solved [40]. 


V. INVERSE & HYBRID PROBLEMS OF FULLY 3-D ROTOR-FLOW 


For these problems three approaches have been developed in China. 

1) Method of Mean-Stream Surface. 

It was originally suggested by Wu in 1952 [1] by a Taylor-serise expansion of flow 
parameters in the azimuthal direction as an extension of MSLM [ 3 0> 2]. It was improved; 
numerically elaborated and applied to design by Zhao et al. in Refs. [41; 42]; where 
an annular constraint condition is set up> which must be satisfied to ensure that the 
hub/casing walls are axisymmetric. 


2) Universal Image-Space Approach. 

In Ref. [43] Liu developed a universal image-space theory of hybrid problems for 
fully 3-D potential flow; which is a generalization of the image-plane approach of 
Ref. [13], Applying tensor calculus and Stokes theorem; the basic flow equations are trans- 
formed into the following integro-dif ferential equation system for the stream functions 
TJTj and l ^e image space CFIg. 14b>* 


° r , . ,) dOVJQ jMUj) 

tfPjV.VJ h ld(iU 2 3 ), ?(?,$'), d (V,¥) 


( 15 ) 


(16) 


and Eq.(5); where V 1 are the contravariant components of the velocity / in a body-fitted 
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nonor thogonal curvilinear coordinate system ^ l (Fig.l6a). Similarly. a corresponding 
potential function formulation by integro-d i f ferent ial equations of this theory has been 
also presented by Liu et al. [44]. 

A similar method for solving 3-D hybrid problems was put forth by Chen et al. [45]. 
using, however, a second-order PDE formulation: 



This equation was solved numerically by the method AF2. the multigrid technique and 
the artificial density in the transonic region. To greatly simplify the numerical 
solution^ the inverse problem is modified in such a way that the physical contravariant 
velocity components ^(assuming that coordinate is the streamlike line) rather than 
the fully velocity yf is prescribed on the blade surface. An axial compressor rotor was 
redesigned and improved by this method as shown in Figs. 15 & 16. 

3) Unified Variable-Domain Variational Approach 

The variable-domain variational approach[21. 22] has been extended by Liu to hybrid 
problems for fully 3~D incompressible [46]. compressible [47] and transoni c [48] flows in 
rotors. Only one of the VPs is given below for reference: 

j 2 (m;;o= i fffi {/- Mw- f njf i V + l . 

(V) 

where the boundary integral term L takes different form for different problem type. Note 
that the variable-domain variation of J 2 should be taken at the unknown boundaries A^* 
(blade surface)and A 2c f (free trailing vortex sheet). Corresponding numerical solutions to 
incompressib le(H^xD)-prob lem of Mizuki's centrifugal compressor [5 1] and to compressible 
( [Hc+D]XD)-prob lem of an axial turbine stator have been obtained by a novel FE with self- 
adjusting nodes in Refs. [49] and [50] respectively and are shown partly in Figs. 17 & 18. 

VI. OPTIMIZATION OF AERODYNAMIC DESIGN OF BLADING. 

1) Optimization of S 2 -Flow. 

A basic and very important problem in this context is the optimal flow typeCi.e. 
optimal radial distribution of flow) in bladings. This problem was first studided by a 
variational method by Liu [52] and later by Xue[53] and Lu[54]. Recently, this problem was 
treated by an optimal control method by Gu * Miao[55]> so that various inequality design 
constraints can be accounted for. 

2) Optimization of 2~D Cascades. 

Theory of optimization of cascade profile shape can be founded on the basis of one 
of the following flow models. 

i) Simplified modeKLeFol l. Citavy) 


Suction side of 
a single airfoil 





Boundary layer 



7 ! 


Inviscid external flow 
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ii) Refined model ( Li u [57-59] ) 


o 1 


Cascade 



Inviscid core flow 

1 

L_ 

n 

r 

Boundary layer 


Perhaps. a third. most accurate flow model based on Navier-Stokes equations incorporating 
some turbulence models should be also tried. 

Two possible mathematical modelsC formulations) can be used here, namely: a) mathema- 
tical programing problem, b) optical control problem. 

2-1. Optimization based on the simplified ftow model. 

The problem of determining the optimal velocity distribution along the suction side 
was formulated as an optimal control problem with some inequality design constraints 
Ce.g. separation-free, maximum or minimum velocity limit, etc) and solved by a heuristic 
grapho-analy t ical method by Liu [56. 63]. This approach is followed by Wang [60] in the 
design of an axial ventilator. Having calculated the optimal velocity distribution along 
suction side and specified a reasonable airfoil thickness distribution, a Hc-problem was 
solved by the image-plane method given in Refs. [5.6] to yield the optimal airfoil shape. 
A similar method, with some modifications, for optimizing 2~D compressor cascade was 
presented by Hua & Chen [ 6 1 ] > where a method for estimating the airfoil circulation was 
given and the optimal airfoil shape was obtained by MSLM[20. 30]. 

The optimal velocity distribution along the suction side on a general streamsheet of 
revolution was obtained by Zou[62]. using Mager's transformation of turbulent boundary 
layer. 

In Ref. [63] some generalizations of the LeFolt's optimization theory of blades were 
given by Liu & Wu to accommodate different objective functionals with more general 
constraints. 

2~2. Optimization based on refined model. 

In Refs. [57. 59] Liu suggested a new theory of optimal 2-D cascades based on the 
above-mentioned refined flow model, in which this problem has been formulated as an 
optimal control problem with multiple inequality design constraints on controt-and phase- 
spaces. Two typical optimal control problems were considered: cascade with minimal losses 
and cascade with maximal loading (circulation), and a duality theorm between them has 
been proved theoret ical ly. so that it is sufficient to study only the cascade with minimal 
losses. This theory has been generalized to 2-D compressible flow and to a 3-D axial-flow 
rotor by Liu [58] . The essential feature of this theory lies in its capability of handling 
a wide variety of practical design constraintsC from stress, vibrational, cooling and te- 
chnological considerations) in a unified manner so as to make the optimal solution surely 
feasible and suitable for use in practice. 

2-3. Local optimization of transonic cascades 

Jiang et al. suggested a numerical method for weakening shocks in transonic cascades 
by local optimization of airfoil shape [64]. The airfoil contour segment near the shock is 
represented by a cubic parabola with free coefficients a^.a^.aj.a^. Then the Mach number 
just before the shock Ms is minimized with respect to a*. 

3) Optimal Design of Diffusers. 

The optimal design of 2-D diffusers was considered by Gu & Ji[66] using optimal 
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control for searching optimal wall shape that Maximizes the pressure recovery of diffuser. 
A more general optimization problem of diffuser was put forth by Liu et al. in Ref. [65], 
where not only the wall shape but also the wall suction distribution that maximize the 
pressure recovery without boundary layer separation are sought by optimal control method. 

An 1-D optimal design method for turbine annular axial-radial exhaust diffuser was 
presented by Ling & Jin [67] based on an approximate loss model. The pressure recovery 
coefficient of an optimal diffuser designed in this way has been shown higher than the 
conventional one by 7* by expriment. 

4) Other Optimization Problems. 

Making use of the Parson's number and the concept of optimal reaction degree, Yao 
presented a method for optimizing aerothermodynamic parameters in one-and multi-stage 
steam turbine design. Some guidelines for optimal design of long twisted blades are given. 

In Ref. [69], based on the diffusion factor and equivalent diffusion ratio, the 
optimal solidity problem of 2-D compressor cascades is formulated by Liu as nonlinear 
programing problems, whose analytical solution in form of simple formulae is very 
convenient for practical use. 

A simple approximate method for determining the optimal relative azimuthal position 
of two blade rows in tandem cascades is suggested by Wu & Feng [70] using a simple total 
pressure loss model. 

VII. MISCELLANEOUS INVERSE, HYBRID * OPTIMIZATION PROBLEMS IN ENGINEERING SCIENCES 

In Ref. [71] the finite element method is generalized by Liu * Zhao via variable- 
domain variations in such a way that the nodes are movable. It allows both optimal grid 
and flow field to be cogenrated simultaneously and naturally using directly the VPs of 
aerodynamic problems. 

The inverse problem of heat conduction with unknown boundary was handled by Liu & 
Zhang [72] using Ritz's and FEM based on Variable-domain VPs. An alternative method for 
solving this problem was suggested by L i u [ 7 3 ] by introducing an image plane JY* (T and 
Y’stand for tranrformed temperature and heat stream function respectively). An interest- 
ing invariance property of the nonlinear inverse problem solution with respect to varia- 
ble conductivity is pointed out. An example is solved by FEM based on a pair of comple- 
mentary extremum principles. 

The inverse and hybrid problems of free surface flow under gravity over a dam are 
posed and handled by Liu via VPs in an image-plane [74] and VPs with variable domain 
in the physical p lane [75, 76] . 

In Ref. [77] Liu suggested a novel problem in gear theory — optimal meshing(i.e. 
optimal tooth profile) of spur gears and its variational theory. An analytical solution 
to the optimal meshing with minimal friction losses has been obtained and it has been 
revealed that the cycloidal gearing with radial tooth profile on the lower half tooth 
height used widely in watches and clocks can be regarded approximately as a practical 
gearing with maximal efficiency. 

VIII. CONCLUDING REMARKS 

Research on inverse, hybrid and optimization problems is of great theoretical as 
welt as practical importance in engineering sciences. To our experiences, the three new 
approachesC image-space approach and VP-based approach, especially its variable-domain 
variational variat, for inverse and hybrid problems; optimization approach based on 
optimal control) suggested and intensively developed in China in the last two decades 
have proved to be efficient tools for inverse design and optimization not only in turbo- 
machinery aerodynamics in particular but also in engineering sciences in general and 
deserve further development and application to practice. Design engineers and industry 
will surely benefit a lot from them, if a complete set of computer codes based on these 


155 

Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences 

fTHDES-Iin Editor: G.S. Dulikravich. Washington D.C.. October 22i2L 1PPL- = 

f— 

approaches can be finished and organized into a coiputerized autonated interactive 
design syste*(something like that of Ref. [81]. 
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Fig. 1 Definition of Symbols for fully 3~D hybrid probleas 
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Fig. 2 Cascade flow and its iaage plane 
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Fi9.ll FE solution of H g -proble« Fig. 13 Airfoil shape by FE solution of 
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Fig. 14 3-D rotor flow and its iaage space 
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Fig. 16 The shape of the blade after and before the 


modification 
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Fig. 18 FE solution of [(H c +D)xD] -problem of an axial turbine stator 
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